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ABSTRACT 

Some transiting extrasolar giant planets have measured radii larger than predicted by the standard 
theory. In this paper, we explore the possibility that an earlier episode of tidal heating can explain 
such radius anomalies and apply the formalism we develop to HD 209458b as an example. We find that 
for strong enough tides the planet's radius can undergo a transient phase of inflation that temporarily 
interrupts canonical, monotonic shrinking due to radiative losses. Importantly, an earlier episode of 
tidal heating can result in a planet with an inflated radius, even though its orbit has nearly circularized. 
Moreover, we confirm that at late times, and under some circumstances, by raising tides on the star 
itself a planet can spiral into its host. We note that a 3x to 10 x solar planet atmospheric opacity 
with no tidal heating is sufficient to explain the observed radius of HD 209458b. However, our model 
demonstrates that with an earlier phase of episodic tidal heating we can fit the observed radius of HD 
209458b even with lower (solar) atmospheric opacities. This work demonstrates that, if a planet is left 
with an appreciable eccentricity after early inward migration and/or dynamical interaction, coupling 
radius and orbit evolution in a consistent fashion that includes tidal heating, stellar irradiation, and 
detailed model atmospheres might offer a generic solution to the inflated radius puzzle for transiting 
extrasolar giant planets such as WASP-12b, TrES-4, and WASP-6b. 
Subject headings: planetary systems — planets and satellites: general 



1. INTRODUCTION 

The most useful experimental data constraining 
evolutionary models of extrasolar planets and their radii 
[Rp) come from measurements of transiting planets^. 
The radius of a transiting planet is inferred from transit 
lightcurve measurements, which in combination with ra- 
dial velocity measurements remove the planet mass Mp - 
inclination angle degeneracy. Much theoretical effort has 
been undertaken t o model, and then to u nderstand, the 
meas ure d radii (iGuillot et al. Il996l : iBurrows et al.l 
20001: iBodenhei mer et al.l 1200 l l: 'Burrows et alj 
2OO3I: [Bod cnheim er et all 120031: iB araffe et al. 200'a 
Burrows et al.l I200l iFortnev fc Hubbard 200l: 



Bara.ffe et al.ll2004l: iChabrier et ah, 2004 : .Laughlin et all 



2005; 
2007 



Baraffe et al.l I200E 
Fortnev et all 12007 



l2006t IBurrows et al 



Marlevetal.1 12007 



Chabrier fc Baraffa TloOTl : iLiu et al.ll2008l : IBaraffe et all 
200l) . 



As IBurrows et al.l (j2007D have emphasized, custom 
fits are preferred for each planet to "project out" the 
effects of age, planet mass, stellar flux, etc. before 
one can conclude whether a measured transit radius 
can or cannot be reproduced by theory. Indeed, the 
transit radius of an extrasolar giant planet (EGP) de- 
pends on many parameters (see IBurrows et all 120071 for 
a sensitivity study). These include the planet mass, 
the stellar irradiation flux fp, the transit radius ef- 
fect teurrows et alll2003l : IBaraffe et al.ll2003( ). the atmo- 
spheric composition, the presence of heavy elements in 
the envelope or in a central core, atmospheric circulation 



that couples the day and the night sides, the planet's 
age, and any effects that could generate an extra power 
source in the interior of the planet. It had once been 
thought such extra power could be due to obliquity tides 
when the planet is in a Cassini state (jWinn fc Holmai] 
|2005| ). but such a possibility has now largely bee n 
ruled out (|Levrard et al.l 120071 : iFabrvckv et"airi2007f) . 
However, heating due to the pe netration and dissi- 
pation at depth of gravit y waves (IGuillot fc ShownianI 
12002; Sho wman fc G uillot 2002), or stro ng tidal effects 
([Bodcnheimer et al.i.2003 : Liu c t al,. 2008i; iJackson et al.l 
l2008bl l^ldl) have not been eliminated. The latter could 
be due the pumping of eccentricity b y an und etected 
companion (jBodenheimer et al.l [20011: iMardfind [2007). 
We note that it is not known whether the tidal 
heat is deposited pr edominantl y in the convectiv e core 
(lOrilvie fc LinI [200l [Good man fc Lackned [200l or in 
the radiative envelope (|Wuil2005ah . We make the default 
assumption of all planet radius evolutionary modelers to 
date that it is deposited in the core, but the reader is 
encouraged to keep an open mind. 

The initial conditions assumed for our simulations can 
affect the structures of EGPs at young ages (< a few 
10 Myr). As a consequence, a model of planet forma- 
tion would be useful to determine the pl anet's initial en- 
tropy and, therefore , its early radius (Marlev et al.ll2007l : 
iFortnev et al.ll20(j8l ). However, since the ages of the ma- 
ture transiting EGPs discovered so far are a few Gyrs, 
precise initial entropies are not crucial for our study. 
On the other hand, as we show in this paper, if tidal 



effects in the planet or in the star are taken into ac- 
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(jFord et all 120031 : iChatterieeet al.l 120081: iFord fc Rasiol 
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I2008t lJuric fc Tremaind I2008D. or t he operation of 
the Kozai mechanism (iWu fc Murray! 120031: Iwiil 120031 : 
iWu et al.|[2007l: iNagasawa et al.ll2008D do matter. 

There are two classes of discrepancies between 
measured transit radii and theoretical predictions. 
The first is the subset comprised of those EGPs 
larger than default radius prediction s and includes 
HP 209458b (jCharbonneau et al.1 l2G00^ iKnutsoneLalJ 
[200l. TrE S-4 (Mandushcv ct al. 2007! ), WA S P-12b 
(iHebiD et al.| [2009ll . WASP- 4b dW ilson ct al. [20081 : 
iGillon et al.ll2009bD . WASP -6b (iGiilOT i ct al. 2009a), XO- 
3b dJohns-KruU ct al. 2008; Winn ct al. 2008), and HAT- 
P-lb ijBakos et al...2007; Winn et al.,i,2007. : Johnson et all 
|2008[ ). The other class consists of EGPs that are m 
fact smaller than default predictions without den se cores 
and, therefore, seeni to re quire such cores (Guill ot et al.l 



120061 [Burrows et al.ll2007l) . The correlation Guill ot et all 
(|200fi ) and lBurrows et aLl (|2007f ) discovered between the 
inferred core mass and the stellar metallicity, which is 
not used in the modeling, supports both the notion that 
dense cores are present and the core accretion mechanism 
for giant plan et formation. 
Previously, iBurrows et al.l (|2007t l studied the effects of 

an extra heat source (i^tide) in the planet's interior and 
estimated the power ne cessary to explain some of the 
measured radii. iLiu et al. (2008 ) took this a step further 
by incorporating core tidal heating to explain the ob- 
served radii of TrES-4, X0-3b, and HAT-P-lb, but kept 
their or bital paramete r s const ant at currently observed 
values. iJackson et al.l (|2008bl l3ldl) studied the coupled 
evolution of e, a and -Etide of close-in EGPs, but main- 
tained a constant planetary r adius. 

H owever, a, nd importantly, IJackson et al.l (|2008bl l3ldl) 
and lGu et al.l |2003) highlighted the possible role of tidal 
orbital evolution in inflating the radius of an exoplanet. 
In this paper, we verify their original insight by consis- 
tently and simultaneously coupling the evolution of Rp, 

e, a, -Etido, and Fp (-^insolation) in the context of a sophis- 
ticated atmosphere and irradiation model for the planet. 
We take into account the tides raised on the planet and 
the tides raised on the star. We find that if the tides are 
strong enough a planet's radius can undergo a transient 
phase of inflation that temporarily interrupts shrinking, 
and does so at epochs consistent with the few Gyr ages 
of transiting planets. Extremely strong tides fade at a 
very early stage and have a negligible impact on the ra- 
dius at Gyr ages. It is noteworthy that the behavior is 
nonlinear and depends sensitively on the initial orbital 
conditions and on tidal heating parameters. In addition, 
we find that increasing the planet's atmospheric opacity 
accelerates and increases the magnitude of the transient 
effect. If the tides raised on the star are negligible, the 
orbit reaches a final circular equilibrium state. Other- 
wise, there is no final equilibrium, an d the planet might 
eventually plunge into its host star (jRasio et al.l Il996l : 
ILevrard et al.ll200"9l ) . We find that due to an earlier phase 
of tidal heating a planet whose orbit has circularized can 
still have an inflated radius. As an example and proof of 
principle, we apply our formalism to HD 209458b. We 
find a set of tidal parameters, Q'p and Q^, and initial 
orbital parameters, e,; and a,;, that lead to the measured 
values of R 



a, e, even for atmospheric opacities for so- 
lar equilibrium abundances. We find that HD 209458b's 



radius can also be fit with no tides and a 3x to 10 x solar 
planet atmospheric opaci ty at the new and shorter age 
of ~3.1 Gyrs suggested bv lTorres et al.l ()2008l ). However, 
unless the atmospheric abundances can be measured, it 
will be difhcult to distinguish the different predictions, 
particularly given the current ambiguity in Q' . Neverthe- 
less, if there are independent (theoretical?) constraints 
on Q' , and/or precision measurements of the atmospheric 
spectra, one might be able to discriminate between the 
two different explanations. 

In we present our formalism, model assumptions, 
and computational techniques. Section [3] is the central 
section of the paper in which we describe the generic cou- 
pled evolution of Rp, e, a, -EtideZ-Einsoiation when tides 
play a significant role. We describe in §3.1! the typical 
scenario, which exhibits transient radius inflation due to 
tidal heating of the planet. The case of extremely strong 
tides is discussed in ^3.2\ In ^3.31 we address the effect 
of atmospheric opacity and in §3.41 we analyze the in- 
fluence of the tides raised on the star. The latter can 
cause the planet to spiral in. In 21 we apply the model 
to HD 209458b. Finally, in ij5| we summarize our re- 
sults, review caveats, and discuss how such transient tidal 
heating might explain other large-radius EGPs, such as 
WASP-12b, TrES-4, and WASP-6b. 

2. MODEL: FORMALISM, ASSUMPTIONS, AND 
COMPUTATIONAL TECHNIQUES 

The evolution of a planet's eccentricity and semi-major 
axis due to tidal effects depends sensitively on it radius. 
In turn, the evolution of its radius depends sensitively on 
the tidal heating power in its interior associated with the 
corresponding evolution of its orbital parameters. More- 
over, a planet's radius evolution is sensitive to the degree 
of stellar irradiation, which is directly tied to the planet- 
star distance. Therefore, to perform simultaneous orbital 
and radius evolutions, all the relevant equations must be 
coupled and atmospheric boundary conditions that vary 
systematically with changing irradiation regimes must be 
incorporated. We have established the tools necessary to 
self-consistently accomplish such calculations and in this 
section we describe our methods. 

We assume that the planet has a spherical gaseous 
H2,He envelope, an d use the equation of state of 
ISaumon et al.l ()1995l ). The helium mass fraction (Y) is 
set equal to 0.25. The effect of reasonable variations in 
Y on the radius of an EGP is small. Except in the atmo- 
sphere, we assume the planet is fully convective and that 
its envelope contains no heavy elements. Though our 
formalism allows it, in this paper we ignore the posssible 
eff ects of an inn e r den se cor e. As shown, f or inst ance, 
in iGuillot et aD (|2006f ) and IBurrows et all (|2007f l. the 
presence of heavy elements in either the envelope or the 
core decreases the total radius of the EGP. To model the 
evoluti on of the planet , we use the Henyey evolutionary 
code of IBurrows et al.l (jl993l . [l997b . with boundary con- 
ditions that incorporate realistic irradiated planetary at- 
mospheres. Using COOLTLUST Y, a variant of the spec - 
tral atmosphere code TLUSTY (|Hubenv fc Landll995( l. 
we precalculate grids of atmospheres for various values of 
Teff and surface gravity g. These provide the associated 
entropy S in the convective region of the planet. Invert- 
ing the relation to obtain Tg/ / {S, g) yields the energy flux 
escapmg the interior fcrTj^^^" : IBurrows et al.ll2003f l. We 



precalculate these grids not only for a given gravity and 
Teff, but also for each orbital pair (e,a). The stellar 
flux at the planet is taken to be the time-averaged mean 
during an orbit. Thus, we have atmospheric and evolu- 
tionary boundary conditions which follow the evolution 
of the planetary orbit. The stellar spectrum is inter- 
polated at the actual effective temperature and gravity 
of the star fro m the Kurucz stellar atmosphere models 
()KurucJl994D . Note that we assume the stellar spectrum 
is constant during evolution. We include the "transit ra- 
dius effect," which accounts for th e fact that the tran- 
sit radius is an im pact parameter (|Burrows 
IBaraffe et al1[2003l) . 

In our calculations, we assume that tidal heating oc- 
curs entirely in the convective interior of the planet 
and that the evolutionary process starts a few Myr af- 
ter the star's formation. Therefore, we suppose that 
the p rotoplanetary disk has dissipated (|Goldreich fc Saril 
|2003[) and that any ensuing chaotic coU isional pe- 
riod of planet-planet scattering has ended fFord et al.' 
2003- Juric & Tremainc 2008; Chattcricc et al. 2008; 
Ford fc Rasioll2008l : INagasawa et al.ll2008l) . Importantly, 



I2008clldl ial: iFerraz-Mello et all 120081 : iBarnes et alll2009f) : 



we presume that, after these early formation, migration, 
and dynamical phases, an interesting subset of close-in 
EGPs are left with high values of (> 0.2) and small val- 
ues of the "initial" semi-major axes, (<0.1 — 0.15 AU). 
{ai is still larger than the "final" values currently ob- 
served.) These assumptions are necessary for tidal ef- 
fects to be of interest as a possible explanation for the 
large planetary rad ii observed in a subset of cases and are 
not unreasonable (iFord fc Rasio 2008; Nagasaw a et al.l 
l2008f ). We neglect stellar and planetary obliquities, as- 
sume that the planet's spin is synchronized (is tidally 
locked) with its orbital period, and that the star's spin 
rate is small compared with the orbital mean motion. A 
rough estimate of the synchronization time (rgync) for a 
close-in EGP gives a va lue between one and a few x 10 
Myr (iGuiUot etaI][T99l) . We assume that equilibrium 
tides have a constant lag angle for any frequency and that 
Q' is independent of orbital period. Another approach 
is to assume a constant time lag fo r any fr e quenc y, as 
was done in the pioneering work of iDarwi see 
also Wu^ 119811. iRasio et al.l Il996l . lEggleton et"arni998l . 
and iLevrard et al.l 120091). Higher-order e formulations, 
though developed (iMard Ung fc Lin"2002t lOgilvie fc LinI 
12004 iDobbs-Dixon et a l. 2004; MardlinllioOTD , intro- 
duce numerous other uncertain approximations concern- 
ing tidal processes. Given our current limited knowl- 
edge of tidal dissipation/heating in EGPs, we feel the 
approach we have employed is acceptable. 

With all these assumptions, the equations of tidal evo- 
lution of eccentricity e, semi-major axis a, and tidal 
heating rate i^tido, to se c ond order in eccentricity, are 
jGoldrcich fc Sote ri [l966[: IKaulal [l968l: IPeale fc CassenI 
1978; Murray fc Dcrmotri999^,'Bodenheimer et al."200U 
2003; .Gu et al.. 2004. : Mardling 2007 : Jackson ct all 
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where Kip, and A'2p» are constants defined by 
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G is the gravitational constant, and Mp, M^^^Rp, i?* 
are the masses and radii of the planet and star. Note 
that the planet radius is time-dependent {Rp (i)), but the 
star's radius is assumed to be constant. The evolution 
of e and a are due both to tides raised on the planet 
(due to the star) and tides raised on the star (due to 
the planet). Equations IJ) and ([2]) indicate that e and a 
can only decrease. This is because we assume the stellar 
spin rate is low. This is different from, for example, the 
Earth-Moon system, for which the tides raised on the 
Earth increase the semi-m ajor axis, while the tides r aised 
on the Moon de crease it (|Goldreich fc Sotei) 119660 . As 
iMardlingj (120071 ) has pointed out, eqs. ©, ©, and ^ 
don't include the higher-order terms beyond and these 
terms should be important for high e. However, in this 
paper we focus on the generic character of the planet- 
orbit coupling and its potential role in explaining the 
subset of EGPs with anomalous radii. In this light, and 
given the many remaining ambiguities in tidal heating 
theory, our formalism should be adequate to address the 
central phenomena. 

The Q's in eqs. (P), (O, and ^ are the tidal 
dissipation factors [Q'^ in the planet and Qt in the 
star), given by Q' — 3Q/ 2fc2, where Q is the specific 
tidal dissipation function (Goldreich 1963") and fca is 
the Love num ber (Love 1927; Goldreich fc Soteil 119661 : 
lOgilvie fc Lin|[2007[ ). Smaller Q's result in higher rates 
of tidal dissipation. For Jupiter, we have Q' k, A Q 
(|Gavrilov fc Zharkovlll977D and estimates of Q are ^6 x 
10^-10*^ iG oldreic h fc SotCTlfl96HYoder fc Pealdll981h . 
lOgilvie fc Li n (200J) have provided a theoretical motiva- 
tion for values of Q near ~ 10^. For a synchronized 
short-period EGP, w e have only crude estimates, but 
lOgilvie fc (I2004D prop ose values of - 5 x 10*^ - 10^ 
and lJackson et all (|2008bl lcl) propose values of Q'^ ~ 10^-^ 
and Q'p ^ 10^-^. The latter are based on the statistical 
distribution of initial values of e and a the authors obtain 
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after integrating eqs. ((T]) and ^ backwards in time at 
constant Rp. For terrestrial planets, the physics of tidal 
dissipation is different and it is thought that Q' > 10^ 
(iDickev et al.lll994l; iMardling fc' Lin 2004; iJackson et all 
12008a'). Most of the numerical values for this factor are 
empirically determined. For Jovian planets and EGPs 
there are various ongoing effor ts to improve our under- 
standing of tidal dissipation (lOgilvie fc Liril 120041 : IWul 
I2005airbl : [Goodman fc Lacknej|2008D . but there is as yet 
no comprehensive theory. Moreover, theory has yet to 
determine where in the planet the tidal heat is deposited. 
If it is not deposited in the convective core, but in the 
radiat ive atmosphere, our results might need to be al- 
tered. iGoodman fc LackneH ()2008D have studied dynam- 
ical tides in an isentropic fluid body, focusing on iner- 
tial waves excited by scattering from the rigid core of 
a planet, and conclude that these waves might dissipate 
most of their ener gy in the convective r e gion. This is also 
the suggestion of lOgilvie fc LinI (|2004[ ). IWul (|2005a[ ). on 
the other hand, suggests that a large fraction of the tidal 
heat could be deposited in the outer regions, perhaps in 
the radiative zone. However, all modelers who have ex- 
plored the possible effects of tidal heating on EGP radius 
evolution have assumed that the heat is deposited in the 
convective interior where it can be redistributed almost 
instantly. For specificity, we stick to this convention for 
this study, assume that Q' is a property of each body, 
and assume it is constant during the integration. 

Since radius evolution must account for stellar irra- 
diation, an informative quantity is the ratio between 
the tidal heating rate -Etido and the insolation rate. 



insolation — 



stellar point 
ratio is: 

-Etidc 



Fp, where Fp is the flux at the sub- 
To second order in the eccentricity, this 



olatii 



63 



m: 
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(6) 



where a is the Stefan-Boltzmann constant and is the 
effective temperature of the star. As is clear from cq. 
([nj, this ratio is a very stiff function of a, and a mod- 
erately stiff function of Rp. We integrate forward in 
time equations H]) and ^ , along with the Henyey equa- 
tions of planetary structure and radius evolution which 
incorporate tidal heating given by eq. ([3]). We start at 
various arbitrarily specifled initial eccentricities and 
semi-major axes and also specify Q'p and Qt, keeping 
them constant during an integration. The strong non- 
linear coupling betwe en e and a, already emphasized by 
IJackson et al.l (|2008cl ). is made all the more so in our 
more general formalism by the variation in Rp and the 
stiff dependence on Rp in eqs. ((T]) and 

From eqs. ([T]) and we can estimate the relative 
contributions of the tides raised on the planet (due to the 
star) [{de/dt)p, (da/dt)p] and raised on the star (due to 
the planet) [{de/dt)^, (da/dt)^]. We obtain: 



{de/dt)j: 
(de/dt)^ 

{da/dt)f 
{da/dt)^ 



28 /M, Y 
25 \Mp ) 



7e' 



57 2 
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Mp 



Rp 
i?* 

Rp 



Q'p 
Q'p 



(7) 



(8) 



The factor [7eV(l -I- 57eV4)] in eq. ([5]) ranges from 
~ 0.46 (e = 1) to zero (e = 0). If we assume that 
the EGP has the mass and the radius of Jupiter and 
orbits a sun-like host star, then the ratio given by eq. 
([2]) is (ep/e,) ~ 14(Qt/Qp) and the ratio given by 
eq. dl]) is (dp/a,) ~ 88 (Q'JQ'p) [eV(l 57eV4)] . For 
HD 209458b, the factors 14 and 88 become 73 and 457, 
respectively. Thus, if the star is as dissipative as the 
planet {Q'^ ^ Q'p) or less dissipative {Q'^, > Q'p), then 
the evolution of the eccentricity is due mainly to the ef- 
fect of the tides raised on the planet. However, if the dis- 
sipation inside the star is high enough {Q'^ <C Q'p), then 
the effect of the tides raised on the star has a compara- 
ble or even greater influence. Moreover, once the orbit is 
circular (e = 0), the tidal effect on the planet ceases (eq. 
[3]), but since the planet and star are not synchronized 
the tidal effect on the star does not. As a consequence, 
the planet's orbital decay continues, though on a differ- 
ent timescale and it can eventually inspiral into its host 
star. Incidentally, when e = 0, cq. ([2]) can be integrated 
analytically to yield: 



a = oo 



1 - 



117G1/2 Mp 



757 (^-*o) 



-1 2/13 



(9) 



where oq is the semi-major axis at any time to after 
the orbit has cir cularized. A similar formula is given in 
'Goldreich' (1963). This formula demonstrates that spi- 
ral in accelerates when the star is more dissipative (Q'^ 
lower), has a bigger radius i?,, has a lower mass Af,, or 
if the planet has a higher mass Mp. 

3. THE GENERIC COUPLED EVOLUTION OF THE 
PLANETARY RADIUS, THE ECCENTRICITY, AND THE 
SEMI-MAJOR AXIS 

In this section, we present and analyze the generic re- 
sults of interest that have emerged from our calculations 
of the simultaneous evolution of the radius and orbital 
parameters of transiting (close-in) EGPs when tidal ef- 
fects are included. The models assume various values of 
Q'p and Q'^ , as well as representative initial values of the 
orbital parameters and , and are chosen to highlight 
various possible behaviors that might be germane to the 
explanation of the anomalous radii of a subset of the 
measured transiting EGPs. Moreover, for discussions of 
a "generic transiting system" we have employed the prop- 
erties of HD 209458 and HD 209458b (listed in Table [J). 
However, it is important to point out that the values of 
the relevant parameters are specific to each individual 
planet-host star system. In each of the following sub- 
sections, we illustrate our conclusions with a four-panel 
figure (Figs. [Tl[2l[3l and [4]) in two rows, with Rp(t) and 

e(t) on the first row, and a{t) and (i?tidc/£^insoiation)M on 
the second row. 

3.1. A Baseline Scenario 

For clarity's sake, we here neglect the influence of tides 
raised on the star (see tj3.4p . equivalent to setting Q'^ 
oo. Figure [1] depicts the simultaneous evolution of the 
planet's radius Rp{t), eccentricity e(t), semi-major axis 

a{t), and power ratio (i?tidc/£'insoiation)M and demon- 
strates how tides raised on a close-in planet might induce 
a transient phase of radius inflation, with consequences 
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on Gyr timescales. In this baseline case, Q'^ — 10^'^ and 
a,; — 0.075 AU. We plot curves for different initial eccen- 
tricities (e, = 0.65,0.60,0.55,0.40,0.30,0.20) and in Fig. 
[T] use atmospheric opacities for solar-metallicity equilib- 
rium chemical compositions. For comparison, we include 
two curves depicting cases with no tides and, therefore, 
no tidally-induced migration, zero eccentricity, but two 
different initial/final values of the semi-major axis (black 
dotted: a = Oi = 0.075 AU; black dashed: a 0.047 AU, 
the current position of HD 209458b). 

Typical behavior is illustrated by the curves on Fig. 
[TJ for which d = 0.65, 0.60, 0.55. The very first episode 
of radius shrinkage is followed by a transient period of 
expansion, after which the radius resumes shrinking at 
a progressively slower rate. Let us describe the case 
for which Ci = 0.65. When = 0.65, radius expan- 
sion starts at the age of ~0.13 Gyr (Rp ~ 1.48 Rj), 
almost in phase with the increase in the power ratio 
(~ 3 X 10~'^). The peak of the radius is reached at 
~0.85 Gyr (Rp ~ 1.78 Rj), and roughly coincides with 
the peak of the power ratio (~ 8 x 10~^), reached at 
^0.80 Gyr. The corresponding time lag of ~50 Myr is 
not universal, but is generally shorter than the charac- 
teristic thermal cooling timescale. The latter depends 
more directly on diffusion through the thick atmosphere, 
while the former is more dependent on core tidal heating, 
to which the radius responds more directly. The average 
rate of radius expansion from the starting point of this 
phase to the peak is roughly 0.43 Rj Gyr~^. The cor- 
responding average rate of increase in the power ratio is 
^ 7 Gyr^^. After the peak in the power ratio is achieved, 
it is followed by a drop at a rate of ~4 x 10^^ Gyr^^. The 
rate at which Rp decreases is comparable to the rate at 
which it increases, until an age of ^-^1.2 Gyr, after which 
it progressively flattens and tends to zero, with an aver- 
age value between 5 and 6 Gyr of ~0.01 Rj Gyr~^. It is 
important to reemphasize that all the variables e, a, 

and -E'tidcZ-E'insoiation are interdependent and evolve con- 
sistently by mutual influence. Thus, as suggested by eq. 
([6|) , the peak in the power ratio is the result of the com- 
bined evolution of e, a, and Rp, which in turn depends 
on this ratio. During the radius inflation phase, e and a 
decrease at a rate that increases in absolute value. The 
peak value of the planet's radius corresponds closely to 
an inflection point in the evolution of e and a. After this 
inflection point, the absolute value of the rates of change 
of both e and a decrease, finally tending to zero. 

While the circularizing orbit is gradually moving closer 
to the star, the stellar irradiation flux Fp is increasing 
and its role in stanching heat loss from the planet and 
slowi ng radius shrinkage is strengthened (|Burrows et al.l 
I2OOOD . Once orbital equilibrium state has been achieved, 
tidal effects disappear and Rp cont inues to evolve due t o 
radiative losses from the surface (jBurrows "erall[200l . 
Therefore, due to an earlier episode of tidal heating a 
planet's orbit can currently be circular, and yet its radius 
can be larger than the default evolutionary theory would 
predict. In other words, a zero eccentricity orbit and, 
therefore, the absence of current tides, does not preclude 
an inflated radius due to the earlier action of tides. We 
will see in §3.41 that this equilibrium orbital state might 
not persist if we incorporate tides raised on the star and 
the corresponding Q'^ is small enough. Note that even 



if the maximum ratio £'tidc/-E'insoiation can appear quite 
small (e.g., ^ 10~^), it can have a lar ge effect on th e 
orbital parameters and radius evolution ()Liu et al.ll2008D . 

For lower values of e^, the duration of the transient 
phase increases, the maximum Rp achieved decreases, 
and the epoch of peak Rp shifts to older ages. The peak 
disappears altogether for values of Cj below between 0.40 
and 0.30. However, this does not mean that tidal effects 
are no longer important. For lower values of e^, tidal ef- 
fects still slow the decrease in Rp. All else being equal, 
in particular for a given a^, the lower the initial eccen- 
tricity ei , the slower the evolution of the eccentricity to 
zero and the slower the evolution of the semi-major axis. 
In other words, the lower the value of Ci, the slower the 
evolution of both e and a and the higher the final a, for 
a given a,;. As for its impact on radius evolution, a lower 
Ci implies a weaker, but longer lasting, tidal influence. If 
ei is small enough, the planet no longer receives enough 
tidal heat to inflate. All these outcomes are straightfor- 
wardly understood as consequences of eqs. ([T]), ([2]), ([3]), 
and ([6]), with Q'^ 00. Thus, a lower value of implies 
lower values of £^tidc and, therefore, less power to inflate 
the planet or reduce its rate of shrinkage. Clearly, it also 
implies a lower initial power ratio iStidcZ-Einsoiation and a 
slower rate of initial decrease of e and a. Furthermore, 
numerical integration of the set of evolutionary equations 
shows that this initial trend for e and a persists at later 
ages, even if the radii and power ratios "invert." By this 
we mean that the radius at around 1 Gyr when = 0.65, 
for example, can be much larger than the radius when 
Ci — 0.20. However, after ~3.5 Gyr it then decreases 
much faster and ends up smaller. Nevertheless, there is 
a small, but long lasting, tidal heating effect even in the 
case of Ci — 0.20. Eventually, the final radii are larger 
if the final semi-major axes are smaller , consistent with 
the conclusions of Burrows et al.l ()2007[ ) — for a circular 
orbit, the closer the planet is to the star, the higher the 
insolation flux, and the larger its radius at a given epoch. 

Comparing the case without tides, but with a ^ ai ~ 
0.075 AU (black dotted curve on Fig. [1]), with those 
with non-zero values of e,; and the same , demonstrates 
that the planet's radius would always be larger when 
the planet can be tidally heated. However, even without 
tides the radius of a planet already at its current position 
(e.g., a = 0.047 AU, for HD 209458b) can be larger than 
the radius of a planet which experiences tidal heating. 
Fig. [1] demonstrates that, if the planet starts at a larger 
orbital distance, tidal heating does not necessarily and 
universally result in a larger planetary radius at a given 
age. 

3.2. Extremely Strong Tidal Effects 

As described in ^3.11 tidal heating can have a sig- 
niflcant, at times non-monotonic, impact on a close-in 
EGP's radius evolution. This is a key conclusion. How- 
ever, contrary to what one might naively expect, ex- 
tremely strong tidal dissipation in a planet can have a 
negligible effect on its late-time radius. Figure [2] depicts 
the simultaneous evolution of a planet's radius Rp(t), 
eccentricity e{t), semi- major axis a(t), and power ra- 
tio (£'tidc/^'insoiation)[i] undcr such circumstances. For 
comparison, the same two tide-free cases depicted in 
Fig. [1] are also plotted in Fig. O The only differ- 
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ence with what was discussed in < j3.1l is that we assume 
Q'p = 10^ °, instead of = 10^-^ We still assume 
Qi — 0.075 AU and ignore stellar tides. Models with 
only three of the initial eccentricities from Figure [1] are 
shown (ci = 0.40,0.30,0.20), sufhcient to demonstrate 
our point. Note that the plots focus on an earlier inter- 
val (2 Gyr). 

In this case, after only ~10 Myr, whatever the ini- 
tial eccentricity, Rp inexorably shrinks and the associated 
curves are much closer. This is despite the fact that in 
the early stages of evolution the power ratio, for = 0.40 
for example, is two orders of magnitudes higher than in 
the baseline scenario of §3.11 — 10~^ compared with 10"'^. 
However, Fig. [2] indicates that such large additional core 
power has had only a small effect on the radius of the 
planet at the typical ages of observed transiting EGPs, 
i.e. at a few Gyr. This is because the orbit circular- 
izes very rapidly and the tidal power ratio disappears 
extremely quickly. The final semi-major axis is reached 
in less than ^^0.1 Gyr and the eccentricity is damped in 
less than ~0.8 Gyr. If we look at ages earlier than ~1 
Gyr (not shown), we do see a transient phase of radius 
inflation for the case = 0.40. The peak occurs very 
early (~10 Myr), is of very short duration, and though 
the planet's radius can reach ~3 Rj, this transient phase 
lasts less than ~20 Myr. The radii, when there are tides 
(e.g., for a = 0.40,0.30,0.20), are close to one another 
simply because the final positions are close. 

However, we suggest that the rapid onset of significant 
heating, followed by the muting of a significant effect on 
Rp at later times, might be an artifact of our initial con- 
ditions — such a combination of low Q'^, high e^, and 
small Ui may not have been allowed to establish itself 
without a more gradual feedback on the evolution that 
we, curiously, witness as an impulsive response. Before 
Ci would have been allowed to get that high or before 
would have been allowed to get that small, tidal effects if 
Q'p = 10^ would no doubt have come into play. This con- 
clusion does, however, depend upon the processes, and 
their timescales, that result in high and in our initial 
Ui. If early migration and eccentricity pumping can oc- 
cur on very short timescales, an early, almost impulsive, 
tidal response is possible. 

3.3. Coupling Enhanced Atmospheric Opacity with 
Tidal Heating 

If one increases the atmospheric opacity, due either to 
an increased metallicity or to the possible effects of pho- 
tolysis and/or non-equilibrium chemistry, the planet re- 
tains heat better and can maintain a large r radius longer. 
This effect was previously explored by IBurrows et al.l 
()2007( ). who concluded that an enhanced opacity delays 
radius shrinkage and can lead to larger EGP radii at a 
given age. They did their study for fixed circular orbits 
and pres ented resu l ts for solar and 10 x solar atmospheric 
opacity. iLiu et ahl f2008) also examined this effect, but 
at solar, Sxsolar, and lOxsolar opacities and explored 
the potential effects of tidal heating (at fixed, non-zero 
eccentricity). However, both studies kept the semi-major 
axis fixed. Here, we investigate the effect of enhanced 
opacity, but include orbit evolution as in §3.11 and §3.21 
For specificity, we focus on atmospheric opacities equiv- 
alent to those for chemical equilibrium abundances at 



solar and 3 x solar metallicity. 

Figure [3] depicts the simultaneous evolution of the 
planet's radius Rp(t), eccentricity e(i), semi-major axis 

a{t), and power ratio (£'tide/£'insoiation)M for both solar 
(solid) and 3 x solar (dashed) atmospheric opacity. We 
present our results for an illustrative subset of eccen- 
tricities [ci = 0.60,0.55,0.40), ignore tides raised on the 
star, and use the same parameters, Q'^ = 10^'^ and Oi = 
0.075 AU, employed in §3.11 The black curves are refer- 
ence models with no tides and a ~ ai ~ 0.075 AU. The 
black, dashed curve for 3 x solar with no ti dal effects reca- 
pitula tes the results already obtained bv IBurrows et al.l 
(|2007( ). At the age of '^l Gyr, the radius increase effect 
is 5%. 

The case for = 0.60 and with tidal effects is par- 
ticularly illustrative ^. With larger atmospheric opac- 
ity, the peaks in the radius and power ratio during the 
transient phase increase in magnitude, narrow, and shift 
towards younger ages. Note that when tidal effects are 
included in the comparison of the consequences of en- 
hanced opacity, the results are not simply monotonic. 
For Ci = 0.60, an enhanced opacity leads first to a larger 
radius for ages up to ~0.9 Gyr, then to a smaller ra- 
dius for ages up to ~2.8 Gyr, after which the planetary 
radius is again larger. Importantly, the higher the at- 
mospheric opacity, the faster the eccentricity evolves to 
zero and the faster the semi-major axis tends to its final 
value. This is a consequence of the larger values of Rp 
possible with enhanced opacity and the stiff dependence 
of -Etide on Rp. Interestingly, the final value of a is the 
same as that reached for the lower atmospheric opacity. 
Why do we end up with the same final equilibrium circu- 
lar orbit, as indicated in Fig. [31 whereas in i j3.1l different 
values of resulted in different final orbits? This behav- 
ior can be derived directly from eqs. H]), ([2]), and 
([6]) . For a given planet / star pair and the same values of 
Q^, e^, and a^, the amount of orbital energy transferred 
to the planet's interior is the same. The increase in the 
atmospheric opacity decreases the rate of heat loss and 
the planet's radius remains larger, longer. But this in- 
crease does not involve any transfer of energy inside the 
planet's interior. This is why for identical initial orbital 
energies, the final orbital energies and, therefore, the fi- 
nal semi-major axes are the same. Angular momentum 
conservation can also be invoked to help explain this re- 
sult. Hence, when tidal effects subside, the final radius 
with enhanced atmospheric opacity is indeed larger than 
the one obtained with solar opacity, but the same orbital 
parameters, all else being equal. However, as Fig. [3] indi- 
cates, the transient increase in Rp when the atmospheric 
opacity is larger can be quite dramatic. 

3.4. The Particular Role of Tides Raised on the Star 

Now, we examine what might happen when we include 
tidal dissipation in the star and set Q'^ equal to vari- 
ous non-trivial values. Figure |4] depicts the simultaneous 
evolution of Rp{t), e{t), a{t), and (£'tido/-Einsoiation)M for 
Q'p — 10^-^, Oi — 0.075 AU, and solar- metallicty opaci- 
ties, but this time for three different values of Q'^: 10^ 
(dashed), lO"^'^ (dotted), and oo (solid). We replot from 

^ As Fig. [3] demonstrates, the qualitative behavior for different 
values of is the same. 
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Fig. [T] the three cases with no stellar tides, but with 
e^ = 0.60,0.55,0.40, Q'^ = lO^-^, and a ^ a, = 0.075 AU. 
The solid black curve is for the case with no tides at all 
and a = ai = 0.075 AU. 

Let us focus our discussion on the three cases with 
Ci = 0.60, but for different values of Q'^. Including stel- 
lar tides, the lower the value of Q',, the higher, narrower, 
and earlier are the peaks of the radius and power ra- 
tio during the early transient phase. This behavior is 
similar to that encountered in i i3.3l with enhanced opac- 
ity or in i j3.1l with higher initial eccentricity. This time, 
however, it is caused by the additional terms involving 
Q'^ in eqs. ^ and ([2]). As seen in Fig. [H e and a de- 
crease faster, but at a combined rate such that the ratios, 

eVa"/^ in ^^tidc/^^insoiation (eq. E]) and e^/a^^/'^ in E'tidc 
(eq. [3]), are bigger than in the case without stellar tides. 
Therefore, more orbital energy is dissipated faster in the 
planet's interior. As noted in ij3.1l and ij3.31 the evolu- 
tion is quite non-monotonic and nonlinear. Moreover, 
even when there is no pronounced peak, as for = 0.40, 
though the evolution is smoother, the general trends and 
behavior are similar. 

The case of a = 0.60 and Q'^ = 10^-^ is particularly 
interesting. Figure [4] shows that even after the orbit has 
already circularized (e = 0), the semi-major axis, instead 
of stabilizing at a constant value, continues to decrease. 
In fact, the planet starts to spiral inward at an accelerat- 
ing rate, is eventually tidally disrupted, and collides with 
the star. Equation ([9]) in |j2]gives the analytical evolution 
of a after circularization, ignoring tidal disruption. For 
= 10^ '^ and e^ = 0.60, we intentionally stopped the 
evolution at a = 0.02 AU, which is at roughly the Roche 

1/3 

limit for this system: ~2.5(/9,/pp) i?,, where and 
Pp are the average densities of the star and the planet, re- 
spectively. While the planet is spiraling in. Fig. [4] shows 
that Rp starts to shoot up. 

This late-time behavior after the semi-axis has 
achieved small values and the orbital eccentric- 
ity is ~zero has b e en se e n before by iRasio et al.l 



(ll996D.lLevrard et all (|2009D . I Jackson et all ([2009D . and 
iMiller et al.l (|2009D . However, whether such a dramatic 
effect would obtain depends sensitively on the value of 
Q'^ (Fig. S]) and on whether tidal effects in the planet 
can synergistically force a to achieve such low values that 
tidal dissipation in the star, which does not depend upon 
the eccentricity, can take over on stellar evolutionary 
timescales. It would seem unlikely that the current sam- 
ple of transiting and close-in EGPs are those for which 
we are just now catching a last glimpse before they are 
eaten by the star. The implication for the parent pop- 
ulation of exoplanets would appear extreme. However, 
at this stage, we cannot eliminate the outside possibil- 
ity that, for a subset of very close-in EGPs, stellar tides 
might be important on timescales short compared with 
their corresponding stellar ages. 

4. APPLICATION TO HD 209458B 

The previous sections were devoted to investigating the 
generic behavior of Rp and orbital parameters when tidal 
dissipation in either the planet or the star is at work. Our 
main result is the emergence for a range of values of 
and Oi of a transient phase of radius inflation which tem- 
porarily interrupts radius shrinkage. This phenomenon 



might explain the larger-than-otherwise-expected plane- 
tary radii of an interesting subset of the close-in tran- 
siting EGP population, either as a vestige of an epoch 
of earlier tidal h eating, resetting the evolutionary clock 
()Liu et al.l l2008f ). or as a current episode of significant 
tidal heating. An objective is to find a set of realistic 
values of Q'p, Q'^, e^, and for which the currently mea- 
sured values of Rp, e, and a can be simultaneously ex- 
plained. Other factors, such as the atmospheric opacity, 
the possible presence of a dense core, and the effect of 
a large heavy-element burden in the planetary envelope, 
also come into play, and custom fits for each EGP are 
necessary. 

In this section, as an example we perform such an ex- 
ercise for HD 209458b. However, the general procedures 
and the generic features apply to all large-radius EGPs 
(e.g., TrES-4, WASP-12b, and WASP-6b) for which tran- 
sient tidal effects such as we have found in this paper are 
promising solutions. The observational data we use for 
HD 209458b are given in Table [1] Figure [5] depicts solar- 
opacity model evolutions of HD 2094578b's radius Rp{t), 
eccentricity e(t), semi-major axis a{t), and power ratio 

(£'tidc/£'insoiation)W- Wc supcrposc On Fig. [5] the mea- 
sured values of Rp, e, and a, with their error boxes. Also 
included on Fig. [5l but only in the top left panel, are 
four curves as reference cases with no tides. The dotted 
curve is for a = Ui = 0.085 AU. The other three non-solid 
curves are for solar, 3 x solar, and 10 x solar atmospheric 
opacities and for a = 0.047 AU, the measured value of 
HD 209458b's semi-major axis. 

To fit the observations of HD 209458b, we have ex- 
plored the effects of varying Q'^ and Q'^ from 10^ to 10^, 
from 0.10 AU to 0.055 AU, and e, from 0.80 to 0.20. 
As a red line on Fig. [5l we identify the model with the 
best-fitting parameter set, obtained by trial and error. 
The associated parameters are: Q'p = 10^'^^, Q'^, = 10^ °, 
Qi = 0.085 AU, and e,; = 0.77. To illustrate the ex- 
treme sensitivity to model parameters, we plot curves 
for Bi ranging from 0.79 to 0.72, in steps of 0.01. Thus, 
by invoking this set of initial and planetary parameters 
and evolving the coupled suite of equations for Rp, e, 
and a, we can simultaneously explain the measured ra- 
dius, eccentricity, and semi-major axis of HD 209458b. 
In i j3.11 we described the sensitivity of the evolution to 
Ci. Indeed, when Ci decreases, the transient peaks of 
the radius and power ratio widen, their height dimin- 
ishes, and the epoch of transient heating shifts to older 
ages. The semi-major axis decreases more slowly and set- 
tles at higher values. With Q'^ — 10^ °, the tides raised 
on the star are not strong enough to lead to inspiral of 
HD 209458b into its host star on a timescale less than 6 
Gyr. However, for Q'^ = 10^ '' there is still a slight de- 
crease of the semi-major axis during the 14-Gyr duration 
of the simulations (not shown). This demonstrates that, 
for this choice of parameters, the orbit of HD 209458b 
is effectively stable. For comparison, we have verified 
that values of Q'^ of 10^ ° and 10^-^ lead to inspiral of 
HD 209458b. However, the fit to the measurements is 
not as good. This example demonstrates how sensitive 
the results are to the values of the chosen parameters. 
Clearly, stronger independent constraints on these tidal 
Q's are desirable. 

Since the paper of iBurrows et all ()2007l ). iTorres et al.1 
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(|2008f) have reappraised the age of HD 209458b, from 
b.btlitoS.lto^ Gyr. As a consequence, as Fig.[5]shows, 
the curves for 3 x solar and 10 x solar, but without tidal 
heating, now intercept the error box. We can in principle 
find a better fit if we choose an opacity between 3 and 
10 X solar, but still without invoking tidal effects. Hav- 
ing said that, the introduction of tidal effects enables us 
to fit HD 209458b's radius easily even with solar-opacity 
atmospheres. In addition, we can now in principle pro- 
vide a reasonable explanation for its observed eccentric- 
ity (and/or upper limit) and semi-major axis. More- 
over, if one needed to invoke high-metallicity (not just 
high-opacity) atmospheres to explain HD 209458b's ra- 
dius and, thereby high-metallicity envelopes, a transient 
tidal phase with significant heating now allows even the 
shrinking effect of heavy elements in the envelope to be 
compensated for. This system is a good example of a 
planet whose orbit has almost circularized, but that still 
has an inflated radius. Hence, as we have suggested in 
this subsection and in ^3.1[ the radius of HD 209458b 
may be due to the former action of tides during an ear- 
lier epoch of non-circular evolution. This explanation 
might also be germane to the TrES-4, WASP-12b, and 
WASP-6b planets, as well as to other large-radius EGPs. 

There exists an interesting, if counterintuitive, system- 
atic behavior in models for the simultaneous evolution of 
a, e, and Rp. For a given final position (a, e), an initially 
higher d does not necessarily result in a larger radius. 
The radius evolution is a strongly non-linear function of 
the initial conditions and as Fig. [T] indicates, different 
e^s result in different radius evolutionary curves that in- 
tersect. Comparing radii for models with different initial 
values of a.i and requires one to specify the age of the 
planet. Among the 7 curves on Fig. [5] (for 6^=0.72 to 
0.79), only those for = 0.77 and = 0.76 lie within the 
a, e, age, and Rp error boxes. Moreover, and perhaps 
counterintuitively, at the age of HD 209458b the radius 
for ei=0.77 is smaller than the radius for 6^=0.76. 

Hence, a model with a higher initial eccentricity peaks 
at a higher tidal heating rate, but at earlier ages. At the 
later times at which models with lower initial eccentricity 
peak, they have higher heating rates than the former at 
that same time. In the same vein, for the same e^, start- 
ing with lower Ui results in a higher tidal heating rate at 
earlier ages than the tidal heating rate that eventually 
obtains when starting with a larger a^. 

5. CONCLUSIONS AND DISCUSSION 

In this paper, we have found that if an EGP after early 
migration and dynamical evolution is left in a tight or- 
bit {ui < 0.2 AU) with a modest to high eccentricity 
(>0.2), its subsequent evolution due to tidal dissipation 
can qualitatively alter our interpretation of its measured 
radius. Using a formalism in which the planet's radius 
and orbit are consistently and simultaneously evolved, 
we have found that a transient phase of rapid tidal heat- 
ing and radius expansion can help explain the large radii 
measured for some transiting EGPs even after this phase 
has subsided. This explanation is straighforward for HD 
209458b, but it might also be a factor in the larg e radii 
obseryed for, for example, TrES-4 ([M andushcv ~et al.l 
[2007t iTorrcs et al.' 2008; Sozzctti ct al. 2009), WASP- 
12b {Hebb ct al 2009), WASP-6b LGillon et al. 2009aj), 
and WASP-4b (|WiIson et al.l [20081 : iGillon et al.l l2009bl : 



iWinn et al.l[2009l) . 

We parameterized our models using a range of planet 
and star tidal dissipation factors, a range of initial values 
of the eccentricity and semi-major axis, and three real- 
izations of the planet's atmospheric opacity. Our main 
conclusions are: 

• A giant planet's radius can undergo a transient 
phase of inflation due to tides that temporarily in- 
terrupts its shrinkage and resets its evolutionary 
clock. The upshot is that, for suitable parameters 
(e.g., Q'p and of 105 °-10^-°, a,; less than ~0.2 
AU, and > 0.2), Rp is measurably larger than it 
otherwise would have been, even after tidal heat- 
ing has subsided. We have demonstrated that a 
planet whose orbit has circularized can still have 
an inflated radius due to the former action of tides. 

• Extremely strong tidal heating in a planet will fade 
early in its life. Under such circumstances, the ef- 
fect on Rp will be negligible at the typical ages of 
the observed transiting planets. 

• Higher atmospheric opacities can enhance and ac- 
celerate the transient phase of radius inflation and 
accelerate orbital evolution. 

• The tides raised on the star also enhance and accel- 
erate the transient phase of radius inflation. How- 
ever, even after the orbit circularizes due to tidal 
dissipation in the planet, the orbit does not sta- 
bilize — for small enough values of Q'^, the planet 
can plunge inward, and is then tidally disrupted 
and consumed by the star. 

• Radius and orbit evolution are strongly non-linear 
and stiff functions of the parameters Q'p, Q'^, 
ei, ai, and atmospheric opacity. Custom fits to 
each planet/star system, rather than pre-calculated 
look-up tables, are to be preferred. 

• The parameters (Qp, Q'*fii, 
ai)=(106■5^ 10^", 0.77, 0.085 AU), with a so- 
lar atmospheric opacity, provide a satisfactory fit 
to the measured radius of HD 209458b, while also 
being consistent with its current orbit. In general, 
the higher the atmospheric opacity the less tidal 
heating needs to be invoked to fit its observed Rp. 

• Given its new, younger age, the measured radius of 
HD 209458b can now also be fit without tides, but 
only if the atmospheric opacity is similar to those 
with equilibrium chemical abundances for metal- 
licities greater than ~3xsolar. However, unless 
the atmospheric abundances can be measured, it 
will be difficult to distinguish the different predic- 
tions, particularly given the current ambiguity in 
Q' . Nevertheless, if there are independent (theo- 
retical?) constraints on Q\ and/or precision mea- 
surements of the atmospheric spectra, one might 
be able to discriminate between the two different 
explanations (vestigal tidal effects and higher at- 
mospheric opacity). 

There are numerous caveats that should be borne in 
mind when considering our results: 

1. Our atmospheric boundary conditions assume 
that day side and night side cooling are the 
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same. Credible general cir culation mo dels (GCMs) 
jGoodman 2008; Sho wman erall 2008b.a,d: 
Dobbs-Dixon fc Lin i2008l: iLangton fc Laughlip 
20081: ICho et all l2008t iMenou fc Rausch cr 200j 



might address this issue most usefully. 

2. We have fixed the properties of the host star during 
our integrations. On timescales of a few Gyr, the 
evolution of the star might be relevant. 

3. The tidal dissipation factors Q' are very poorly con- 
strained and are free parameters of our models. A 
firmer theoretical handle on th eir values and depen- 
dences woul d be very josefuyOgilvic &Jin||^Ol; 
IWu ,2005ak iGoodman fc Lackncr 2008) . 

4. We have assumed that the planet's spin is synchro- 
nized (tidally-locked) and have ignored the spin an- 
gular momentum of both the planet and the star. 
This is generally a good approximation, but there 
may be exceptional cases. 

5. We have assumed that the tidal heating power is 
deposited solely in the convective zone, in which 
it is rapidly mixed. If much of the heat is in- 
stead deposited in the radiative atmosphere, as 
might be expecte d if circularization i s dominated 
by Hough modes (jOgilvie fc Lin|[2"007f ). our results 
would need alteration. 

6. We have neglected the obliquities of both the star 
and the planet. 

7. We start our calculations after early disk evolution 
and dispersal, and assume that the anomalously 
large planets, for which we might invoke tidal ef- 
fects, are left after these early phases with high 



eccentricities and small semi-major axes. These 
assumptions are at present almost wholly uncon- 
strained. We may, however, be able to turn the 
question around and someday use the current ob- 
servations to constrain an earlier episode of vigor- 
ous tidal heating. 



We have demonstrated that, given the right param- 
eters and initial conditions, early tidal heating can 
result in large EGP radii, even while leaving the 
planet's current orbital parameters at modest values 
for which tidal heating is not expected to be im- 
portant. These vestigal effects of an earlier tran- 
sient phase might echo into the present to help ex- 
plain the anomalously large radii of a small, but in- 
teresting, subset of transiting extrasolar giant planets. 
Our formalism can easily be applied to oth er inflated 
trans i ting EGPs, such as TrES-4 (Mandu shev et aO 
[20071 iTorres et al.' '2008'; 'Sozzetti et al. 2009), WASP- 
12b dgcbb et al. 2009), WASP-6b (Gillon et al. 2009a 



and WASP-4b (IWilson et all 120081 : IGillon et all 12009 , 
I Winn et al.ll2009| ). and we plan to do so in the near fu 
ture. 
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TABLE 1 

Observational Data of the HD 209458 System. 
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Note. — Data are from lKnutson et al.l ||2007|) . [Torres et al.l (|2008|'| . and IMadhusudhan fc WinnI ||2009|') . 
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Fig. 1. — The Baseline scenario. Shown are the evolution of the planet radius Rp (Rj) (top left), eccentricity e (top right), semi-major axis a (AU) (bottom left), and power ratio 
^tide/^insolation (bottom right) vcrsus age (in Gyr) for a representative planet (taken to be HD 209458b) at solar atmospheric opacity. Q'p is set equal to 10®'^ and is set equal to 
0.075 AU. We neglect tides raised on the star (Q'^ — > oo). assumes values of 0.2, 0.3, 0.4, 0.55, 0.60, and 0.65 (different colors). Also shown in the left panels are radius evolution curves 
with no tides — a = (dotted) and a = amcasurod = 0.047 AU (dashed). For e; = (0.65,0.60,0.55), the evolution undergoes a transient phase of radius inflation. This phase does not 
appear for lower values of (0.40,0.30,0.20), but the radius increase effect at later times is still in evidence. The orbits eventually achieve circular equilibrium states — the higher the e^, 
the closer to the star. See the discussion in i|3.1l 
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Fig. 2. — The same as in Fig. [T] but for the case of very strong planetary tidal effects (at a low value of Qp = 10^) and for the first 2 Gyrs. Also shown on the two left panels are 
two reference radius evolution curves incorporating no tides — one for a = (dotted) and one for a = Qnioasurcd = 0.047 AU (dashed). Models for only a subset of values of (0.20, 
0.30, 0.40) are depicted. The tidal effects fade at very early stages, so that their effect on radius evolution is quite small at the observed ages of transiting planets, a few Gyr. The orbits 
circularize in less than ~0.8 Gyr. See the text in ^3.2l for a discussion. 
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Fig. 3. — Same as in Fig. [T] but depicting the effect of atmospheric opacity. Models for opacities with equilibrium chemical abundances at solar (solid) and Sxsolar (dashed) metallicity 
are compared. The initial eccentricities are 0.40, 0.55, and 0.60 and they are distinguished by different colors. The planetary tidal dissipation factor (Q'p) is set equal to 10®'® and 
the tides raised on the star are neglected (Q^ — > oo). a; is set equal to 0.075 AU. Also shown on the top left panel are two radius evolution curves with no tides — for a = a, at solar 
(solid) and Sxsolar (dashed) equilibrium opacities. The final orbital state for the two opacities is the same, but that with higher opacity reaches it earlier. See the text in i|3.3l for a more 
complete discussion. 
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Fig. 4. — Same as Fig. ^ but portraying the effect of tides raised on the star. Models for values of of oo (solid), lO'^ O (dashed), and lO^-^ (dotted) are compared, is set equal to 
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Fig. 5. — Example fit of the observed values of Rp, e, and a for HD 209458b. The measured Rp, e, and a are shown with error boxes. The best fitting curve is for = 0.77, Q'p = 10° °"', 
and Q', = 10^". This model is shown in red. We have used a,; = 0.085 AU. To demonstrate the strong dependence on e^, we include models for Ei between 0.72 and 0.79. Also shown in 
the top left panel are four radius evolution curves (in non-solid black) that ignore the effects of tides. They are for a = a^, but also for a = ameasured = 0.047 AU at solar, Sxsolar, and 
lOxsolar atmospheric opacities. The best fitting curves without tidal effects are for >3xsolar. See the text in [|4]for a discussion. 



